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The fitting formula for the location of the first acoustic peak in the matter power spectrum is 
revised. We discuss the physics that leads to baryon acoustic oscillations: the recombination history, 
the tight coupling approximation and the velocity overshoot effect. A new fitting formula is pro- 
posed, which is in accordance within 5% with numerical results for a suitable range of cosmological 
parameters, whereas previous results yield deviations of up to (15 — 20)%. The crucial improvement 
turns out to be the accuracy of the recombination history. 
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I. INTRODUCTION 

Before recombination the Universe is dense and highly 
ionized, and baryons and photons are tightly coupled by 
Thomson scattering. The pressure of the Cosmic Mi- 
crowave Background (CMB) photons opposes gravita- 
tional collapse and leads to acoustic oscillations. In fact, 
during this phase, the amplitude of perturbations in the 
baryon density cannot grow, but they oscillate with a 
slowly decaying amplitude. After recombination, baryons 
decouple from radiation and the oscillations are 'frozen 
in'. Because baryons represent a significant fraction 
of matter, cosmological perturbation theory [1, predicts 
that these acoustic oscillations of the baryons (BAO's) 
are imprinted on the late-time matter power spectrum, 
leaving features analogous to the acoustic peaks in the 
CMB power spectrum. The BAO's have indeed been ob- 
served in the large scale galaxy distribution [2]-[4], and 
they are one of the main observational goals of recent 
and upcoming surveys. 

Numerical calculations of the recombination history 
are available thanks to, e.g., the RiCO [5] or the Rec- 
FAST [5] codes for the recombination history. The lat- 
ter, in particular, is the default model used in the Camb 
Boltzmann-solver code [7]. It does reproduce the results 
shown in p and is a fast approximation to the detailed 
calculations described in [8 with some updates discussed 
in [2| and with the Compton coupling treatment of [TU] . 
The numerical calculation of the recombination history 
is much more time expensive than employing analytical 
approximations or fitting formulae, like the ones used 
by [llj which are an improvement of the fits presented in 
Appendix D of [1]. Nevertheless, an accurate computa- 
tion of the recombination history turns out to be a sig- 
nificant step for the evaluation of the location of troughs 
and peaks in the transfer function. 

At early times, before recombination, baryons and pho- 
tons behave as a single 'tightly coupled' fluid because 
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Thomson scattering, which couples electrons and pho- 
tons, is much more rapid than the expansion of the uni- 
verse ts ^ H~^, where tg = (cTTne)"^ is the photon 
scattering time scale (i.e. the mean time between Thom- 
son scatterings) and the expansion time scale of 
the universe. Here ctt is the Thomson cross section and 
rie is the number density of free electrons. Since scat- 
tering is rapid compared with the travel time across a 
wavelength, we can expand the perturbation equations 
in powers of the Thomson mean free path Xs = ts = k^^ 
over a wavelength A oc k~^, i.e. k/k, where k = aaxne 
is the differential optical depth. To the lowest order we 
obtain the tight coupling approximation (TCA) [T^. A 
more rigorous definition and treatment of the TCA can 
be found in [13j , while [Uj analyzes the second-order ap- 
proximation in the inverse Thomson opacity expansion. 
In [15 formulas for the TCA to second-order are derived 
independently and tested. These results are implemented 
by the Class Boltzmann code. 

It can be demonstrated that at large scales the trans- 
fer function is governed by density perturbations, which 
oscillate roughly like a cos (fcs) where s is the sound hori- 
zon at decoupling, see e.g. [T2] and the Appendix D of 
[1]. However, the corresponding velocity perturbation 
dominates in the small scale limit. When the oscillations 
are released at decoupling, baryons move kinematically 
according to their velocity and generate a new density 
perturbation '11]. This 'velocity overshoot' effect is re- 
sponsible for the fact that the transfer function, for suf- 
ficiently large ks, actually behaves like sin(fcs). 

In this paper we derive a fitting formula for the loca- 
tion of the peaks and troughs in the matter power spec- 
trum by matching the solutions for the matter density 
perturbation before and after decoupling. We obtain a 
form which is consistent with the one proposed in |11] for 
the position of the first peak. However, our fit is tested 
considering recent cosmological parameters and it uses 
an improved recombination history. The latter turns out 
to be an important amelioration which lets us achieve a 
significantly better accuracy than the previous fit. Even 
though one can compute the positions of these peaks nu- 
merically with the help e.g. of CAMB-code, we believe 
that an analytical fit has its merits as it helps us to see 
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immediately what effect a variation of cosmological pa- 
rameters will have and since it gives us a better under- 
standing of the physics involved. 

The paper is organized as follows. In Section [ll] we 
review the physics leading to BAO: the ionization history, 
the TCA and the velocity overshoot effect. In Section |TTT] 
we give the fitting formula and we compare the results 
with pjj. We state our conclusions in Section IV Some 
details of the derivations are moved to an Appendix A. 



II. BARYON ACOUSTIC OSCILLATIONS 

We review the physics leading to BAO, in order to 
highlight the main results that allow us to locate the 
peaks and troughs in the transfer function. 

We use the notation of [T7]. In particular, t denotes 
the cosmic time and rj conformal time such that adrj — 
dt, where a is the scale factor. An over-dot indicates 
a derivative w.r.t. the conformal time (') = d/dr]{). 
We also use of the notation R = 3pf,/4p^, where the 
subscripts b and 7 label the energy density of baryons and 
photons, respectively. Our reference model is a ACDM 
Universe. 



A. Ionization history 

In [S] a calculation of the recombination of H, He I, 
and He II in the early Universe is developed which is im- 
plemented in the publicly available code Recfast. The 
methodology is to calculate recombination with as few 
approximations as possible. One of the main improve- 
ments with respect to previous calculations is that it 
takes into account that the population of excited atomic 
level depart from an equilibrium distribution. Indeed, 
recombination is not an equilibrium process. Simplified 
analytical calculations or approximate fitting formulae 
for the recombination history are too crude to give good 
approximations for the location of peaks and troughs in 
the matter transfer function as we will discuss in §III[ 

Recfast models hydrogen and helium atoms as effec- 
tive three-level systems, including some artificial correc- 
tions via conveniently chosen fudge factors ^ . An alter- 
native approach is the RiCO code 5J, which smoothly 
interpolates the ionization fraction on a set of pre- 
computed recombination histories for different cosmolo- 
gies. The inclusion of previously neglected physics leads 
to changes in the ionization fraction at the (2 — 3% level 
in some redshift regions. However, throughout this pa- 
per, we run the Camb code using Recfast version 1.5, 
which gives a sufficiently accurate recombination history 
for our purposes. 

On the other hand, we treat the late-time reionization 
of hydrogen and helium via the fitting formulae proposed 
in the appendix of [T8]. 



B. Tight Coupling Approximation 

Here we derive an analytical solution for the baryon 
density perturbation in the tight coupling approximation 
(TCA) in first order perturbation theory using the WKB 
approximation valid for a slowly varying R, inside the 
sound horizon at decoupling given by s ~ ch'^^ rjdec- Here 
c^P^'^ is the sound speed of the baryon-photons fluid (de- 
fined in appendix [a] Eq. ( A5 )), and rjdec is the decoupling 
time. Hence, we must keep in mind that this approxima- 
tion is valid only for sufficiently large k. For the range of 
our interest this is fine since acoustic oscillations concern 
relatively small scales, of the order of lOO/i"^ Mpc 0. 

We perform our calculation in the uniform curvature 
gauge, the differences between the variables calculated 
in different gauges is small on sub- horizon scales |19j . 
Furthermore, all the physical observables must be indeed 
gauge invariant. So, in terms of the density perturbation 
in the uniform curvature gauge, Dg \YI\, the general tight 
coupling solution for the baryon density perturbation is 
given by (see appendix |A] for a derivation) 
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(3/4) Dg™^ is determined by the adiabatic initial 
condition and ^ {k, rj) is the Bardeen potential [19] . We 
have introduced the (comoving) sound horizon rg (77) = 
d(c]'' {(), i.e., the distance that a wave can travel in a 
time rj. During the tight coupling phase the baryon den- 
sity perturbations undergoes harmonic motion following 
roughly a cosine mode with an amplitude that decays in 
time as (1 -h i? (r?))"^^''. 

To show that this solution follows a cosine mode, a 
simple analytical approximation of the Bardeen potential 
^ (k, ri) can be obtained by writing the Bardeen equa- 
tion in the case of adiabatic perturbations for a mix- 
ture of perfect fluids (photons, baryons and CDM). On 
super-horizon and sub-horizon scales one finds [19j . re- 
spectively, 

*x«i(fc,f?) = *o(fc) , (2a) 
*.»i(fc,'7) = -3*o(fc)^^ , (2b) 

where the initial metric perturbation \l/o (fc) is constant 
in time and x = k cf^^ drj. To derive these relations, 

we also assume ~ (^ci''^^^ — 1/3. Here Cs — P/p, 
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where P and p denote the total pressure and energy den- 
sity of the baryon-photon fluid, respectively. The latter 
approximation means that, since the WKB approxima- 
tion requires slowly varying R, we suppose i? ~ over 
an oscillation period. This implies that we are also ap- 
proximating the equality epoch as roughly the decoupling 
epoch. 



Computing the integral in Eq. ( lb I we obtain an an- 



alytical approximation for the baryon density perturba- 
tion in the tight coupling limit. Neglecting the small 
contribution from a; < 1 in the integrand of Equation 
(lb), using Equation (2b I and i? ~ 0, wc obtain 
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where Si(^) = Jq dxsinx/x ^nd Ci ) = 
— dxcosx/x the sine and cosine integral 

functions. Finally, we consider modes x ^ n and use 
approximately the asymptotic form Si (x) n/2 and 
Ci (x) — valid for a; ^ 1, we obtain 

I{x) ~ 7r~^ sin(a;) . (4) 

We summarize our result for the tight coupling approx- 
imation in the form 
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Making use of the perturbed Einstein constraint equa- 
tions and of the Friedmann equations to rewrite ^'o , we 
obtain 
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Deviations of -Dg*^''^''' from the cosine mode decay like x~'^ . 
In Fig. [T] wc compare the pure cosine mode with the full 
approximate solution given by Equation (|6| . For x > Stt 
the deviation from the cosine mode is negligible, and only 
for a; < 7r/2 the integral term is dominant. 

In conclusion, the peaks of the tight coupling solution 
for baryon density perturbations D^^^'^'^ [k, rj) closely fol- 
low those of cos (fcc^^Tj) for small scales, say cc > 37r ~ 10, 
while the few first peaks may exhibit deviations due to 
the integral term. In particular we expect large devia- 
tions from the cosine mode for a; < 1. This consideration 
allows us to find a formula to fit the position of the peaks 
and of the troughs of the matter power spectrum. 




FIG. 1: The solid line shows Equation p|, with 

^gT^ ~ 1 ^^'^ -^(■^) evaluated as in Eq. ([3b|. The 
dashed line is cos a;. The curves are in good agreement 
for a; > 37r and deviations are significant only for x < n. 



C. Velocity overshoot 



The velocity overshoot effect can be explained by not- 
ing that decoupling is close to equality, rjeq < rjdec- Be- 
fore 77eg ~ Vdem baryons are tightly coupled to pho- 
tons and their velocity is governed by the dynamics of 
photons, which are the dominant component. Indeed, 
Vf, = sin(fcs), where ^ indicates 'oscillates as', 

is the tight coupling solution in the limit k ^ oo (see 
Appendix [A} Eq. (A2)). When rj > rjeq, the energy den- 
sity of photons pj becomes smaller than the matter en- 
ergy density pm- Furthermore, for rj > rjdec baryons are 
no longer coupled to radiation. This implies that for 
V > Veq ~ Vdec, baryons no longer follow the photon 
velocity. As we shall see, the baryon velocity after de- 
coupling, Vb (jj > rjdec) ^ COS {ks) is almost exactly out 
of phase with (77 > r/dec) ^ sin (ks). 

This can be shown by matching the solutions for the 
baryon density perturbation before and after decoupling. 
As derived above, the adiabatic initial conditions for an 
inflationary model select the cosine mode for the baryon 
density perturbation tight c oupl ing solution on small 
scales, which is given in Eq. (0 for kc'p''^r] > n. Let 
us indicate Dgb the solution after decoupling; we match 
it to the tight coupling solution. 



Dgbik,ridec) = D'gY'\k,r]dec) , 
Dgbik,ridec) = D^g^,'''\k,rjdec) ■ 



(7a) 
(7b) 



After decoupling baryons evolve like CDM. The evolution 
of Dgb can then be evaluated by considering the Bardeen 
equation for a mixture of non-interacting radiation and 
matter fluids in a matter dominated epoch that, neglect- 
ing the decaying mode, yields "if = 'i>{k, rjdec) constant 
in time. Using equations (7a) and (7b) as initial condi- 



tions and denoting the present time by 770, on small scales 
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ks ^ 1 we obtain 

-'^^[gik.ridec) + ridec9{k,ridec)] , (8) 



The dominant term here 



[km) 



^oiVo/Veq)^ which comes from the baryons falhng 
into the gravitational potential of dark matter. In 
addition, we have a growing function oscillating like 
A:ssin(fcs) plus a correction due to the ^Po-term which 
slightly affects the period of the oscillations, see 
discussion below Eq. 

To better understand the expression in square brack- 



ets, let us consider Equation ( 5b ) , which for ks > tt shows 
that 

= -18 7r-ifcscos(fcs) , (9) 

This suggests that the position of the troughs and peaks 
in the matter power spectrum may differ slightly from 
those of ks sin (ks) and this difference is proportional 
to which is in turn proportional to Q,mh^ ■ Since 
the pre-factor is small, we can approximate ks sin(fcs) + 
ekscos{ks) by fessin(fcs + e) + O(e^). With this, we ex- 
pect that the positions of the troughs and the peaks in 
the matter power spectrum are approximately given by 



fc„ = ^ (1 + /3„ 



(10) 



where = 3, 7, 11, . . . for the troughs and n — 5, 9, 13, . . . 
for the peaks, and where /3„ is a parameter that takes into 
account the correction which affects mainly the lowest 
knS and which can be fitted by comparing with numerical 
results. Since ^E'o oc —Dg and g + rjdccg < 0, we naively 
expect e > and hence /3„ < 0. We shall see, however, 
that the approximations made for the decoupling redshift 
actually lead to /3„ > 0. 



III. FIT OF THE ACOUSTIC PEAK POSITIONS 



Eq. ( 10 1 allows us to localize the troughs and the peaks 
in the matter power spectrum. We finally want to derive 
an explicit form for the sound horizon at decoupling s, 
defined as the comoving distance that a wave can travel 
prior to decoupling tdec- 



cp'''^ (1 + z) dt' 



(11) 



The sound speed of the photon-baryon plasma is given 
in appendix [Aj Eq. ( A5 1 . 

This integral can be computed exactly if we neglect the 
contribution of dark energy to z, which is a very good ap- 
proximation for the redshifts we are interested in. The 
subscripts 6, c and m refers to baryons, CDM and non- 
relativistic matter (baryons plus CDM), respectively; we 



define the density parameter ujx = 0.x for the species 
X. The subscript 7 refers to photons while the subscript 
r refers to the density in relativistic particles at the time 
of equal matter and radiation, which probably also com- 
prises three types of neutrinos. We consider, Wc, h, cUm 
as independent cosmological parameters, keeping the first 
two fixed and varying the latter. We then write the re- 
maining parameter as LUb = cOm — ^c- This yields 

^ h p dx 

HoV^ Jl+Zd„ {x + r){xLOr + UJm) 

Ah 



(12) 



i?o is the value of the Hubble parameter today and 
r = {\ + z)R — 3wfc/(4a;^) is the r-parameter defined 
in .2Qj- It is worth stressing that here Zdec denotes the 
baryon decoupling redshift. In fact, even if radiation 
decoupling epoch is z* ~ 1100, the scattering rate of 
electrons is sufficient to keep the matter temperature 
r„i equal to radiation temperature T-y down to redshifts 
z ~ 100. Then, to extract BAO's information is nec- 
essary to integrate up to the baryon decoupling epoch, 
also known as drag epoch. However, for z < Equation 




(A5), which gives c^'^'^ during the tight coupling regime. 



no longer holds. Nevertheless, we note that in the Rec- 
FAST code Tm is set to (plus a small correction) until 
a typical switching redshift z ~ 850, after which the full 
evolution equation for Tm are considered [T^. We find 
convenient for the fit to artificially shift the baryon decou- 
pling redshift to the Recfa st ' switch' value, Zdec = 850, 



still considering Equation (A5) for 

the Zd 



Furthermore, 



value is almost independent of the cosmological 
parameters within the range of our interest [21j . 

This approximate Zdec treatment force us to reconsider 



the interpretation of Equation (10). In fact, we extrapo- 



late this formula to correct also the uncertainties related 
to Zdec- As we will see, this implies in particular that 
/?„ > in our fit. 

In |llj a fitting formula for the matter transfer function 
of a CDM plus baryon Universe can be found. The cur- 
vature and also the cosmological constant are neglected. 
Since the latter do not contribute significantly to the 
sound horizon at decoupling, this approximation is still 
valid in a ACDM Universe. In Figure [2] we compare 
the first peak positions ^i.pfc/fcf'p^' evaluated approxima- 
tively as 5tt/2s. The wavenumber ki^pk is calculated by 
using Eq. ([l2| with Zdec in accordance to Recfast, while 



ki'Jk' calculated according to the sound horizon at de- 
coupling employed in Eisenstein and Hu (1998) [TT]. The 
parametric formulas for s lead to an agreement within 1% 
for the first peak position when considering the same re- 
combination redshift as used in [llj . Instead, if we use 
the value consistent with Recfast in Equation (12 1, we 
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FIG. 2: Comparison of first peak position for h = 0.70, 
LUm — 0.13, LUb = 0.02. ki^pk is evaluated by using Eq. 



(12) for s and the z^ec employed in '11 (solid line) and 



in Recfast (dashed), fef^f- is the fit proposed in [TT| . 



find that in [TT] the position is systematically overesti- 
mated (in terms of k) by about 10%. Furthermore, the 
full fitting formula proposed in [TT], accounting also for 
the fim-correction, yields disagreements up to 15% with 
respect to the numerical results obtained with Camb. 



A. Fit of the acoustic peak positions 

Let us discuss, for illustrative purpose, the fit of the 
first three troughs and peaks in the matter power spec- 
trum, FigurejS] As we are not now interested in precision, 
the fit is compared to a numerical code 22 that agrees 
with Camb within about 5%. 

It is clear that the form of the fit is adequate to repro- 
duce the numerical results, but let us consider Table |lj 
which reports the fit parameters obtained, to check our 



expectations. From Equation (10 1 we see that the rela- 



Order 


Troughs Peaks 


1st 

2nd 

3rd 


0.25 0.07 
0.12 0.08 
0.12 0.10 



TABLE I: /3-correction, defined as Pn'^m- We have 
fixed ujm = 0.144. 

five importance of the /3-correction is given by j3ni^m- As 
explained above, this correction is due to the fact that 
the first nodes of the transfer function slightly differ from 
those of sin(/cs), also for ks > tt, because of the velocity 
overshoot effect; together with the approximation for z^ec 
which is too high, hence an s which is too small, leads to 
f3n > 0. Indeed, as shown in Table|T]a correction of about 
25% is obtained for the first trough. This /3-correction is 



FIG. 3: First three troughs and peaks fit. The circles 
and the crosses are the numerical data for troughs and 
peaks respectively. The dashed and solid lines are the 

fits for troughs and peaks respectively. We fixed 
h = 0.70, Wc = 0.114, ujr = 4.17 • 10-^ uj^ = 2.48 • lO'^ 
and Hf^'^ = 2997.9 h-^ Mpc. 



larger than the other cases, for which is about 10%, and 
it is due to the fact that the corresponding ks — 3tt/2 
is the closest to the critical value ks — it. Actually, we 
also note that the corrections for the troughs are larger 
than for the peaks; this is due to the method used to 
extrapolate the trough and peak positions, but here we 
neglect this detail. 



B. The fitting formula for the first pealc 

The first peak position in the matter power spectrum 
is conveniently fitted by 

fci.pfc = ^ (1 + 0.276 . (13a) 

Inserting the wellknown photon density uj^ and Hq/H we 
obtain for s 

s = 19.9 (w™wh)^'/' log [U i^b, uJm)] Mpc , (13b) 



and 



U iuJb,UJrn) 



1.12 J (1 + 28.18 uj,n) + VI + 35.54 ojb 



1 + 1.12 



Note that the units are Mpc, not h ^Mpc. With this, 
the fit (13a) deviates by less than 5% from the numerical 



results of Camb for the range of cosmological parameters, 
around the values reported in [53], 0.70 h < 0.75, 
0.100 < < 0.130 and 0.0125 < Wb < 0.030. 

In Figure |4] we show the location of the first peak ki pk 
as a function of baryon and matter density parameters. 
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FIG. 4: The location of the first peak in Mpc~^ as a 
function of baryon and matter density parameters. 
Lines of constant ki^pk are indicated. 



As the baryon fraction D,f,/fl„i increases, the first peak 
is shifted to smaller scales, since the sound speed and 
with it s decrease. The value of ki^pk also increases with 
flm, due to the larger contribution of the J7m-term in 



in [231 , with an accuracy of about 5% with respect to 
the numerical results of Camb, whereas the fitting pro- 
posed in [TT] disagrees by up to 15% with Camb within 
the range of parameters we explored. The 5% error es- 
timation takes into account not only discrepancies with 
respect to Camb, but also minor contributions that may 
come, e.g., from neglected non- linear effects, or varia- 
tions of 1% order in the ionization history. The latter 
might be improved by using a different recombination 
code like Rico and, as minor contribution, by improving 
treatment of reionization . 

The fit may be improved on the one hand by better 
approximations of the baryon density D{,g and by im- 
proving the recombination history and on the other hand 
by improving the fitting formula (13a|-(13b| itself. This 



may also imply a more accurate discussion of the baryon 



decoupling redshift, see Equation (12), which here has 
been artificially shifted to Zdec = 850, based on the im- 
plementation of Recfast, independently of cosmological 
parameters. However, we expect that most of the indi- 
vidual effects will lead at best to a 1% improvement and 
recalling that on this level e.g. also the value of the scalar 
spectral index affects the peak position, we believe that 
a better fitting formula would be very complicated and 
to obtain significant improvements in the evaluation of 
the first peak position one has to resort numerical calcu- 
lations which might need to go to second order in per- 
turbation theory. 

Our fitting formula is especially useful for a first esti- 
mate of the effects of changing cosmological parameters 
on the positions of the baryon acoustic peaks. 



Equation ( 13a 



IV. CONCLUSIONS 

Matching the tight coupling approximation, Eq. (|6]), 
to the solutions after decoupling, allowed us to develop 
further the approach initiated in [Ij. This yields an ana- 
lytical formula for the location of the peaks and troughs 
in the matter power spectrum, Eq. ( |10[ ). The formula 
has the same form as the one given in |llj. 

Using the same approximation for the recombination 
history as [11], we obtained results compatible with [11] 
within about 1%, even though we consider very differ- 
ent cosmological parameters, ~ 0.7 as compared to 
r^A ~ which was considered in [11]. This shows that 
the acoustic peak positions are not really sensitive to f^A 
but only to ujm, tJb and of course ritotai- This corresponds 
also to the findings of j_20j. However, considering an im- 
proved recombination history, i.e., using Recfast, the 
fit proposed in [IT for the location of the first peak in 
the matter power spectrum no longer holds. 

This leads us to propose an improved fitting formula 
for the position of the first peak obtained by running 
Camb, see equations (13a|, (13b) and Figure |4] The fit 
yields the location of the first peak in a convenient range 
of cosmological parameters around the values reported 
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Appendix A: Tight coupling approximation 



In order to derive Eq. (la), we consider the evolution 



of baryon perturbations during the tight coupling regime. 
We follow [21] and [T^, where the evolution of photon 
perturbations is discussed in detail. Since baryons are 
coupled via Thomson scattering to photons, the evolu- 
tion of baryon perturbations is related to that of pho- 
tons by the differential optical depth k — aaxn-e: where 
ctt denotes the Thomson cross section and Tie the elec- 
tron number density. Indeed, the equations governing 
the baryon perturbations evolution read [19j : 

Dgb = -kVb , (Ala) 

Vb + nVb = k-i' + ^ (v^ - Vb) , (Alb) 

H 

where we use R = |^ and V denotes the veolocity per- 
turbation. We also write the first moments of the Boltz- 
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mann equation for photons 



Dg~j — ' 

1 



K = 2fc* + -kD 
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k [V^ - Vb) 



(Ale) 
(Aid) 



Since CDM does not interact other than gravitationally, 
we do not need to consider its evolution here. 



If we take the Umit /i — > oo in Equation (Aid) we find 

Vb = V^. (A2) 

This zero-order tight couphng solution leads to an im- 
portant consideration: during the tight coupling phase, 
perturbations between baryons and photons are roughly 
adiabatic on all scales due to Thomson scattering. 

Using this zero-order result (in 1/k) back in the l.h.s. 
of Equation (Aid I we find the leading order equation: 



k 



K - 14 = V 2* -f -D, 



1 



'97 



K 



Using this in Eq. (|Alb|) we obtain: 
R 



Vb 



1 



k 

i?^^" " 4(l + i?)^^^ 



(A3) 



fc* . (A4) 



Differentiating Eq. (Ala I and using Eq. (A4| to replace 
Vb we find: 



D 



R 



1 + R 



mvb 



R 



4(1 +i?) 1+R 



We use again Eq. (Ala 



that until photons and 



to substitute kVb and the fact 
^aryons are tightly coupled, the 
adiabaticity condition Dg^ = ^Dgb holds [T9|. Then, we 
have 



D 



R 



9b 



l + R 



-no 



9b 



3 (l + R)' 



-D. 



R 



9b 



1 + R 



Using R oc a, the comoving Hubble parameter writes 
H = a/a = R/R. We write the sound speed of the 
photons plus baryons system as 




1 



v/3(l + i?) 



(A5) 



where we also used Pb = and = p^/3 for the baryon 
and photon pressure, respectively. 

We finally write the equation for the baryon density 
perturbations as 



Dgb + if^Dgb + k' [ci^'^y Dgb = F{k, t) , (A6) 
where we have defined the forcing function 



F{k,t) 



R 



1 + R 



k^'^{k,t) 



(A7) 



Eq. (A6) represents damped, driven oscillations of the 
baryon density perturbation. The second term on the 
left-hand side is the damping of oscillations due to the 
expansion of the universe. The third term on the left- 
hand side is the restoring force due to the pressure. The 
forcing function is governed by the gravitational potential 
perturbations. These oscillations are called 'acoustic os- 
cillations' since, as in acoustic waves, the photon-baryon 
fiuid cannot simply collapse under gravity because of the 
restoring force provided by the pressure which leads to 
oscillations. 



To obtain an analytical solution to Eq. (A6), we first 



find the solutions to the homogeneous equation through 
the WKB approximation |16j . valid for slow varying R, 
inside the sound horizon at decoupling given by s ~ 
cp^'^^Vdec- Then we obtain a particular solution by the 
standard Green's method imposing adiabatic initial con- 
ditions. This yields the general tight coupling solution 
for the baryon density perturbation 



D 



(t.c.) 
9b 



{k,ri) 



D 



{in) 

9b \ \ . 

Eik,r,) 



1 



1/4 



COS {krs 



(A8) 



where 



E{k,v) = {1 + R{v)) 



-1/4 



dc 



ii + R{C)f^ 



^ sin[fc.(,) ...(C)] ^,^(^^^^ 
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